% 一个非常简单的脚本，但是可以用来讲一个有趣的故事
% 假设系统中有两种粒子（以0和1代表）
% 根据统计力学的等概率假设，对于一个孤立系统，所有可能的微观态出现的概率是相同的。
% 然而，那些“有序”的微观态（例如所有的1都在左侧）的个数非常少，因此整体上看，系统总是处于“无序”的混合状态
% Gitee Repo


clc
clear


L = 1;
dx = 0.025;

[x,y] = meshgrid(-L:dx:L);
n = size(x,1);
u = zeros(n,n);


TICK = 1000;
figure

for tick = 1:TICK
    mask = randperm(n*n)(1:floor(n*n/4));

    u(:)=0;
    u(mask)=1;

    clf
    hold on
    axis equal
    surf(x,y,u);

    % 打印单侧粒子1的占比。理论上应该接近50%，即没有显著的聚集
    disp(sum(u(:,floor(1:n/2))(:)/sum(u(:))))

    drawnow
    pause(0.01);
end

